Approximate Schur Complement Preconditioning of the Lowest-Order Nodal Discretizations

Certain classes of nodal methods and mixed-hybrid nite element methods lead to equivalent , robust and accurate discretizations of 2 nd order elliptic PDEs. However, widespread popularity of these discretizations has been hindered by the awkward linear systems which result. The present work overcomes this awkwardness and develops preconditioners which yield solution algorithms for these discretizations with an eeciency comparable to that of the multigrid method for standard discretizations. Our approach exploits the natural partitioning of the linear system obtained by the mixed-hybrid nite element method. By eliminating diierent subsets of unknowns, two Schur complements are obtained with known structure. Replacing key matrices in this structure by lumped approximations we deene three optimal preconditioners. Central to the optimal performance of these preconditioners is their sparsity structure which is compatible with standard nite diierence discretizations and hence treated adequately with only a single multigrid cycle. In this paper we restrict the discussion to the two-dimensional case; these techniques are readily extended to three dimensions.